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Evolution of spatial anisotropy created in the Heavy Ion Collisions at Relativistic Energies (HI- 
CRE) has been studied. The microscopic evolution of the fluctuations has been examined within the 
ambit of Boltzmann Transport Equation (BTE) in an expanding background. The expansion of the 
background composed of Quark Gluon Plasma (QGP) is treated within the framework of relativistic 
hydrodynamics. Spatial anisotropies of the initial state with different geometry have been evolved 
through Boltzmann equation and studied how these anisotropies dissipates during their propagation. 
This is very useful in understanding the presumption that the measured anisotropies in the data 
imitate the initial state effects. Relation between thermal fluctuation and transport coefficients like, 
shear viscosity and thermal conductivity have been presented. We have also estimated the angular 
power spectrum of the anisotropies in the appendix. 

PACS numbers: 25.75.-q,25.75.Dw,24.85.+p 


I. INTRODUCTION 

The main aim of the heavy ion collision experiments at Relativistic Heavy Ion Collider (RHIC) and Large 
Hadron Collider (LHC) is to create a new state of matter called quark gluon plasma (QCP). Such a state 
of matter, i.e. QCP may have existed in the early universe after a few micro-second of the big bang. One 
of the motivation to create and study QCP in the laboratory is to understand the state of the universe 
in the micro-second old era. The temperature fluctuation in the cosmic microwave background radiation 
(CMBR) has provided crucial information about the universe when it was about 300,000 years old. This 
information has led to tremendous support to the big bang model of the cosmology. The polarization of the 
photons resulting from the Thomson scattering at the decoupling surface infected by density fluctuation gets 
reflected in the quadrupole moment of the phase space distribution of the incident photon. The temperature 
fluctuation in the CMBR is introduced as a perturbation in the phase space distribution of photons. The 
evolution of this perturbation is studied by using Boltzmann transport equation (BTE) [1] in gravitational 
field with Thomson scattering in the collision term. The linear polarization resulted from the scattering is 
connected with the quadrupole moment of the photon’s phase space distribution. 

We have adopted a similar approach here to study the spatial anisotropy implemented by introducing 
a deviation, df{x,p,t) in the equilibrium distribution function, /*'°^(p). The bulk in equilibrium evolves 
via relativistic hydrodynamics and the evolution of the fluctuations over the equilibrated system, Sf is 
treated within the ambit of BTE. In contrast to this the propagation of perturbation has been studied using 
hydrodynamics in [2] (see also 0 )- Sf may be used to estimate the fluctuations in various thermodynamical 
quantities as we will see below. In the present work the equilibrated background is assumed to be quark 
gluon plasma (QCP) expected to be produced in heavy ion collision at relativistic energies (HICRE). 

The fluctuations in the thermodynamic quantities (e.g. hot spots created in the initial state of the collisions 
i) can be related to perturbations in the phase space distributions in the hydrodynamic limit. The evolution 
of these fluctuations can be analyzed within the ambit of BTE [S] . Fluctuations in thermodynamic quantities 
have been proposed as signals of the critical end point in the QCD phase transition lilZl- Dissimilar 
fluctuations in partonic and hadronic phases in the net electric charge and baryon number may shed light on 
the QCD phase transition in HICRE [8]. Event by event fluctuations in the ratio of positively to negatively 
charged pions may be used as an indicator of QCD transition [3] as well as for understanding the chemical 
equilibrium in the system formed in HICRE m- Evolution of these fluctuations near the critical end point 
has been studied by using BTE m- Kinetic theory approach has also been adopted to study fluctuations in 
particle and energy densities [12]. In Ref. m it has been argued that the perturbations in hot QCP travel 
longer distance to reach the border of the medium giving rise to the possibility of detectable signatures of 
these perturbations. 

The study of the fluctuations in the space time structure of the fireball driven by the fluctuations in the 
position of the nucleons in the colliding nuclei is an important contemporary issue in HICRE. Fluctuations 
in the space-time structure of the system will infect fluctuations in the thermodynamic quantities. How will 
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these fluctuations evolve with time in a hydrodynamically expanding system and how are they connected 
with the transport coefficients are addressed in this work. 

The paper is organized as follows. In the next section we discuss the evolution of fluctuations in a 
non-expanding background and present a relation between energy fluctuation in Fourier space and viscous 
coefficient. Section III is devoted to discuss the progression of fluctuations in a hydrodynamically expanding 
QGP background. Results are presented in section IV and section V is dedicated to summary and discussions. 


II. EVOLUTION OF FLUCTUATIONS IN A NON-EXPANDING BACKGROUND 

In the following subsections we discuss the connection of (5/, deviation of phase space distribution from 
its equilibrium value with the fluctuations in various thermodynamic variables and its time evolution in a 
non-expanding background within the framework of BTE. The phase space distribution function, f(x,p,t) 
of a system slightly away from equilibrium, at time t, position x, momentum p can be written as 

f{x,p,t) = f^°Hp){l + 'i'{x,p,t)}, (I) 

where f^^\p) is the equilibrium phase space distribution function determined by the local temperature, 
T{t, x) and '^{x,p, t) is the fractional deviation from f°{p). 'if{x,p, t) can be used to estimate the fluctuations 
in various thermodynamic quantities in the system. The evolution of dt is governed by BTE, which in turn 
provides the relation between the dissipative effects and fluctuations in the hydrodynamic limit. 


A. Fluctuations in various hydrodynamic quantities in Fourier space 

The energy momentum tensor, of the system under study can be written as: -|- AT^'^, 

where the equilibrium (ideal) part, is determined by /°(p) and the dissipative part, is determined 

by dt, i.e. 

= J d^P^^-^f{x,P,t) ( 2 ) 

where f{x,p,t) is given by Eq. The of the system in equilibrium can be obtained from Eq. I by 
setting = 0. The metric, in Minkowski space-time is taken as = (—1,1,1,1). We assume that 
the momentum, p can be written as pi = prii, i.e. p = \p\h where h is a unit vector and (Pp = p'^dpdfl, dfl 
being solid angle associated with rii which satisfy f dQuirij = AirdijlS and f dflninjUk = 0. 

It is straightforward to obtain various components of from Eq. The deviation of the components of 
the stress energy tensor, ST^'^ from their ideal values can be expressed in terms of the perturbation, 'k as 
follows: 


<51)? = - Jp'^dpdnef^°\p)'i>{x,p,t), 

<5T? = / p'^dpdnpn^f^°\p)'i>{x,p,t), (3) 

6T^ = J p^dp dil — n^rij (p)^'(x, p, t) 

where e = = —po = \/p^ -I- w? and m is the mass of the particle. It is to be noted that the integral, 

/ d^ppnif^°^p) = 0 since in equilibrium the system is isotopic i.e, all directions are equally probable in its 
rest frame, the phase space average of momentum vector then is zero. 

The ideal part of the stress energy tensor in the hydrodynamic limit is given by, 

T^^ = {p + P)U>^U, + Pg^, (4) 

where = dx^/dr is the four velocity of the fluid and r is the proper time. The relations among various 
components of and the thermodynamic variables are: Tg = —p, Tj = —Tg = 0, and T* = Pd*. 
Therefore, p = j p^dpdQ.ef^°\p), is the average energy density. Other thermodynamic quantities like 
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pressure, number density etc can be estimated in a similar way. Equilibrium distribution is isotropic, 
therefore, integration over dfl will give simply 47r. If the fluid is slightly away from equilibrium with space 
time dependent fluctuations in energy density, pressure, velocity, etc., then the system will evolve toward 
equilibrium through dissipative processes. In such situation the components of the energy momentum tensors 
can be explicitly expressed in terms of the thermodynamic variables as: 

To{xz,t) = -{p + Sp{x„t)}, 

T°ixi,t) = -n = ip + P)v., 

T;{x,,t) = {P + 6P{x,,t)}S; + mx,,t), 

= 0 , 


where Vi is the component of the velocity perturbation. We can choose a frame which is moving with 
velocity close to the velocity of the fluid, so that the fluid velocity measured from this frame is small. From 
Eq. we get E* = By using Eqs. i 0 and [^we get, 

6p{xi,t) = -ST^{x„t), 

6T°{xi,t) 

(p + P) 

STj{xi,t) = 5P{xi,t)5] + T,]{x^,t), 


i{Xi,t) = 


( 6 ) 


with 6P{xi, t) = 6T^{xi,t)/3. The shear stress, E* (xi, t) can be expressed in terms of shear viscous coefficient, 
77 as T,j{xi,t) = —and the thermal conductivity (y) is defined through the relation, 
STj = —x(§^ + [13 ■ useful to express these quantities, i.e. various components of in 

Fourier or fc-space because expansion of these quantities in terms the spherical harmonics (Yim) will enable 
to connect the angular scales set by I in terms of k analogous to the determination of angular scale in 
CMBR [T]. In fc—space these quantities are marked by~as: E^(/ci,t) = —iri{IJ^kj + Ujk} — '^Sjk^IJr) and 

6f°{ki,t) = -X ikif{ki,t) + f{ki,t)^ . 

By using Eqs. ^ and the fluctuations in momentum space can be expressed as: 


6p{k„t) = J p‘^dpdnef^°\p)^{k„p,ni,t), 

Vi{ki,t) = - ] - [p^dpdnpnif^°'>{p)^, 

r ^ 1 - 

^]{kut) = J p'^dpdn^{n,nj - -Sij)f^°'>{p)^, 

5P(k,,t) = \ jp‘^dpdn^f^°\p)'^ 


where 47 is the Fourier transform of 47. Now we take the zenith direction along k and then the angular 
dependence of 4'(fci,p, n^, t) can be expressed in terms of angles between k and n. Depending on the symme¬ 
tries of the problem under consideration 4^ can be expanded in a series of suitably chosen basis functions e.g, 
for axial symmetry in terms of Legendre polynomials and in absence of such symmetry it can be expressed 
in terms of spherical harmonics HU- 

The vector component rii and tensor components {uirij — appearing in the expressions for Vi{ki,t) 

and Td^{ki,t) respectively, can be converted into functions of 6 (angle between k and n), by taking contraction 
with suitable tensors made out of the components of fc. If we contract rii with ki then we get kk-fi = k cos 9 = 
kPi(k ■ n) and by contracting {niUj — ^Sij) with {kikj — ^Stj) we get |^(3(fc • h)^ — 1) = |5(3cos^ 9 — 1) = 
|P 2 (^ ■ n), where Pi s are Legendre polynomials. For axial symmetric distribution of p it helps to connect 
different co-efficients of expansion of 47 in terms of Legendre Polynomials with corresponding scalar quantities 
obtained from Vi{ki,t) and E'^/,t) due to orthogonality relation satisfied by P^s. We define the scalar 
quantities like A, 9 and a as [17] . 


A(fcj,t) 


Sp{ki,t) 

P 


6f§{ki,t) 


P 


( 8 ) 
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is the fluctuation in energy density in Fourier space. Similarly we define energy flux, 

ik^STj{ki, t) 


9{ki, t) = ik^Vj = 


Cp + P) ’ 


and the shear stress as, 


{p + P)a{ki,t) = -{hkj - ^6ij)T,j{ki,t) 


(9) 


( 10 ) 


The quantity, 9{ki,t) = ik^Vj originates from the velocity gradient. 9 and a can be expressed in terms of 
the shear viscous coefficient ( 77 ) and thermal conductivity (y) as follows: 


= -^j^ik^Uj{k,t) 


( 11 ) 


and 


9{k,t) = - - j^{k'^f{ki,t) - ikif{ki,t)Ui{ki,t)). 

p + P 


( 12 ) 


Now Eqs. mill] and can be used to obtain the fluctuations in the energy density, pressure and velocity 
in fc-space as: 


a77 1 /■ JO 

4ttJ Jp^dpef(o)(p) 

5P{ki,t) _ 1 f p‘^dp^f^°^p)^{h,p,ni,t) 


P 47r ^ 

0 Ic C 

9{ki,t) = — / dVL{k.h) 


47r 


f p'^dpp'^/ef^^y{p) 
f p'^dp f^°\p)^{ki,p,ni,t) 
/p 2 dp(e + p2/3e)/(°)(p) 


(13) 


1 


a{ki,t) = -— / dn{{k.n)- -) 


p‘^dpf^°^p)^{k^,p,n^,t) 


dlT 


3' Jp‘^dp{e+p‘^/3e)f^^^p) 


To understand the angular scale determined by the multipole number I (as used in the appendix to find the 
angular correlations) we expand 4/ in terms of Legendre polynomials for an axially symmetric distribution 
as: 


^{k,n,p,t) = + l)'^/i{k,p,t)Pi{k.n), 

1=0 


(14) 


where the factor(—z)*(2Z + 1) is used to simplify the expansion of a plane wave form of 'h. The I is related 
to angular resolution of the anisotropies, i.e. smaller angular scale will require larger I and vice versa. The 
temperature fluctuations may be written in terms of 'I' as AT/T = —[dlnf/dlnp]’^ [17) . For simplicity we 
will consider the massless limit, m = 0 which gives e = p. The energy density is a scalar quantity whereas 
the velocity and the shear tensor are vector and tensor respectively and these aspects are also bound to 
reflect in the corresponding fluctuations. Therefore, the orthogonality of P;’s ensure that the fluctuations in 
scalar, vector and tensor quantities are dictated by the coefficients 4>o, 4>i and 4>2 which are obtained by 
substituting 'I' from Eq. [^in Eq. |13|and performing the angular integration as. 


Aiki,t) = Sp{ki,t)/p = J p'^dppf’'°\p)^o{ki,p,t)/p, 
SP{ki,t)/P = J p^dppf^°^p)'^o{h,p,t)/p 

9{ki,t) = 3ik^6T°/{4p) = Jp'^dppf'^°'>{p)'i>i{k^,p,t)/p 
= ^ J P^dppf^°\p)'S2ih,p,t)/p 


(15) 
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where p = f P^dppf^^^ (p). The above set of equations can be written in a more compact form as: 

F{n,k,t) = J p'^dppf^°\p)'^{k,p,ni,t)/p 

therefore, F has the angular dependence of 'f' and consequently F can be expressed as: 

OO 

F{k,n,t) = ^(-i)'(2/ + l)Fi{k,t)Pi{k.n). 

1=0 


with 


Fi{k,t) = j p'^dppf^°\p)'i/i{k,p,t)/eo 


(16) 


(17) 


(18) 


The fluctuations in terms of T/’s are now given by 

A{k,t) = Fo{k,t), dP(k,t)/P = Fo{k,t), 0(k,t) = S/AkFi{k,t), a(k,t) = l/2F2{k,t). (19) 


Using the relation a{k,t) = —Ar]ik^Uj{k,t)/3{p + P), and writing ik^Uj{k,t) = Q{k,t) we get an important 
relation which connects the fluctuation {F 2 ) with the transport coefficient ( 77 ), 


F2{k,t) = pM kp) = -^^e(fc,t) (20) 

3(p + P) 31 s 


where the thermodynamic relation, p + P = sT has been used. In Eq 20 s is the entropy density and T 
is the average temperature. The p appears as a coefficient of 2"'^ rank tensor involving gradient in the fth 
direction of the jth component of velocity, as a result the I = 2 term appears in the expression for 77 in | 20 | 


B. Fluctuations in Fourier space and transport coefficients in relaxation time approximation 


The temperature fluctuations, 6T{9,(p) in CMBR is generally expanded in Laplace series in terms of 
spherical harmonics, Yim{9,<p). The maximum value of I is determined by the angular resolution of the 
detector which can be connected to the wave number (fc) corresponding to the Fourier transform of the 
spatial anisotropy. Therefore, in analogy with fluctuations in the CMBR the spatial anisotropy is studied 
here in Fourier space. But first we briefly discuss it in coordinate space. 

The BTE, p^dfif = C[f] in absence of external force and in the relaxation time approximation (similar 
approximation were used e.g. in [I8l - l20] l used reduces to 


dt € dx^ Tn 


( 21 ) 


for 'll. In Eq[^ tr is the relaxation time. In this work we assume that the system is close to the local 
equilibrium and the collisions between the particles bring the system back to the equilibrium within a time 
scale Tr. For the present work the relaxation time can be estimated as the inverse of the reaction rate of the 
quarks and gluons using pQCD cross sections and Hard Thermal Loop Approximations m- The solution 
of Eq. 21 for a given initial (at time to) distribution, 'I'i„(i:,p, <o) is [H]: 


'S{x,p,t 


to) = 



:it-to)),p] exp 


jt - to) 

Tr 


( 22 ) 


Knowing 'I' it is straightforward to estimate the fluctuation in energy density from the following expression: 

/p2dpdHe/(°i(p)^'(T,p,t - to) 


/p2(ip(iHe/(°l(p) 

The solution l^is useful to study the time evolution of spatial anisotropy of the matter. 


(23) 
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Now we would like to derive a relation between the fluctuation in energy density and transport coefficients 
in a non-expanding system. To facilitate this we write Eq. 21 for massless particles (as the case may be for 
partonic plasma produced in HICRE) in momentum space: 


, -1 M.T, 

—— \-ik(k.n)w = -. 

ot T 

With the help of Eq. Eq. can be reduced to an equation describing the time evolution of F as 

——h ik{k.n)F = -. 

at T 


This equation has the following solution, 

F(k, n, t) = Fik, n, to) exp 


-(- -I- ikp){t - to) 


(24) 


(25) 


(26) 


where k.fi = p. The value of F{k,h,t) can be obtained from its value at initial time, to- Eq. 26 is a general 
expression for the fluctuations in the sense that all the quantities, e.g. A, 6, ct, discussed above at time t 
can be obtained from this expression if their corresponding initial values are supplied. Expanding F{k, n, t) 
as in Eq. 17 and using the orthogonality relations of Pi{p)s we get. 


Fi{k,t) = -e 


1 _ (t-tp) 


OO . + 1 

Y,{-iY"~'H^s + l)F,ik,to) / dfiPi{fi)Ps{fi)e 

8=0 




Therefore, the energy density fluctuation at time t is given by: 

OO „ + l 


1 /‘“T-L 

A(fc,t) = -e-^ V(-^)*(2s + l)E,(fc,^o) / 

^ s=o "'-i 


Taking terms upto s = 2 in the expression for A(fc, t) we get, 

A(fc,t) = -e-^^(-i)«(2s + l)P,(fc,to) / 

^ s=o "'-i 

Performing the integration over p we get, 

A(fc,t) + 3Pr(fc,to){^°^^^^ ~ sinfc(t-to) 


(27) 


(28) 


(29) 


k{t - to) 


k{t-to) {k{t-to)Y 


} 


“ ^o) , ^cosk{t-to) 3smk{t-to) 

2( , o){ + (fc(t-io))2 {kit-to))^ 


(30) 


This is the fluctuations in energy den sity, fro m whic h the fluctuations in temperature can be obtained as 


6p/p = 4,ST/T for p ^ T^. We use Eqs 
coefficients as: 
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and 30 to obtain the density fluctuation in terms of transport 




k{t - to) 


k s T 


k{t-to) fc2 (f-to)^ 


40 77 0 (A:,to) rSinfc(t — to) 3 cos k{t —to) 3 sin A:(t — to) 11 
3 s T k(t — to) k?{t — to)‘^ fc3(t —to)^ 


(31) 

This equation provides the connection of the energy density in Fourier space with various transport coeffi¬ 
cients e.g. thermal conductivity (x) and viscosity {rf). The fluctuation in energy density in position space 
can be obtained by taking Fourier transformation of Eq. |31|as, 


^{x,t) = 

P 


d?k 

(27r)" 


A(fc, t) exp [ik.x) 


(32) 


If the initial (t = to) energy density fluctuation, gradient of velocity, viscosity to entropy ratio and temper¬ 
ature of the system in equilibrium are known then Eqs. and can be used to get fluctuations at any 
time, t > to- Angular correlation function for these fluctuation has been estimated in the appendix. 
































III. EVOLUTION OF FLUCTUATION IN A HYDRODYNAMICALLY EXPANDING QGP 

BAGKGROUND 
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In the previous section we have considered the evolution of the fluctuations in a non-expanding background. 
However, in a realistic scenario in HICRE the system expands very fast due to its high internal pressure. 
Therefore, in this section we include the effects of the expansion on the spatial anisotropy through the 
solutions of relativistic hydrodynamics. 

The space time variation of quantities such as energy density and flow velocity is governed by relativistic 
hydrodynamics. Therefore, we solve the equation: 




(33) 


with the assumption that net baryon number density (ub) is zero at the central rapidity region, hence we 
need not consider the equation di_i{nBU^) = 0. We also assume boost invariance [53] along the longitudinal 
direction and solve the Eq. numerically for initial condition taken from optical Glauber at the highest 
RHIC energy {^/snn = 200 GeV) for Au-I-Au collision. The hydrodynamic solutions [24] for flow velocity 
and temperature (T = [30p/((77r^)]^/'^) have been used to study the space time evolution of the fluctuations 
in an expanding background. 

The evolution of f{x,p) = fo{x,p)+6f{x,p) with 6f = fo'tp is governed by the equation, p^d^f = {p-u)C[f] 
which reduces to the following for an expanding system in the relaxation time approximation: 


/ 9 p d 
dx 


{p°uo - p-u) \ 

P°tr{x) J 


Sf{x,p) 




fo{x,p) 


The solution of Eq. 34 is given by 


Sf{x,p) = D{t,to) 


Sfin{.P,x- + [ B{x- 

Jto P 


(34) 


(35) 


where 


D{t 2 ,ti) = exp 


f dt'Aip, X- - to), t') 

Jti P 


(36) 


with 


and 


For 


A{p,x,t) 


p^Uo{x) — p- u{x) 
P°tr{x) 


B{x, t) 




fo{x,p) 


fo{x,p) = feq 


I 


The expression for B reduces to: 


(37) 


(38) 


(39) 


B{x, t) 


/e(j(l + feq) ^0 I 


{x)u^ 


Pl^\ 


(40) 


We took the relaxation time as Tj^^{x) = \AasT{x) [5T] (we have taken constant value of Us = 0.2 here), 
/3 = l/r(a:), u^^{x) = ( 7 , 7 u) is the four velocity of the fluid and 7 ( 2 :) = vf{x) = (1 — u(a:)^)“^/^. The 
space-time variation of the temperature and velocity fields are taken from the solution of the relativistic 
hydrodynamic equations. To simulate different types of initial spatial anisotropy we choose. 


5f{j), x,to) = Ao exp [-r(1 -I- a„ cos n-d)] (41) 

We have taken n = 2,3,4, 5,... to simulate different initial anisotropy. For numerical results Aq is set as 
unity. 
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FIG. 1: Evolution of the fluctuation in energy density with r coordinate at different t for a non-expanding QGP 
background. 


IV. RESULTS 

We display the spatial variation of the fluctuation in Fig. at different times for a non-expanding back¬ 
ground. We substitute eq. |^in eq. [^with '1'™ as a Gaussian in space at the initial time and evaluate the 
evolution of the fluctuation. For the sake of illustration we take T = 400 MeV resulting in r ~ 0.98 fm/c. 
The results indicate a rapid dissipation and displacement of the peak of the fluctuations with increase in 
time. Now we present the results for a realistic scenario where the background QGP is expanding hydrody- 
namically as described in the previous section. Taking the value of the temperature dependent relaxation 
time relevant for QCD plasma m the evolution of initial spatial anisotropies introduced through (or 
Sf) have been studied. It is to be noted that due to expansion the temperature decreases and hence the 
relaxation time increases which slows down the dissipation. Therefore, the dissipation of the perturbation 
gets slower with the expansion of the system. In Fig.|^the evolution of the initial elliptic spatial anisotropy 
(upper panel), realized by taking n = 2 in Eq. is depicted. The initial thermalization time is taken as 
To = 0.6 fm/c. The evolution is studied upto r = 4 fm/c when the temperature is ~ 200 MeV. For n = 2 
the anisotropy has an elliptic shape having stronger gradient along x-axis resulting in faster expansion along 
X compared to y axis. Therefore, the propagation of the perturbation generates a pattern similar to the 
one generated in water waves by the impact of a stick on the still water surface. This kind of pattern is 
clearly observed in the lower panel of Fig.[^ It is also important to note that the two oppositely propagat¬ 
ing perturbations are strongly correlated which may have interesting observation effects. Moreover, if the 
perturbation is created near the boundary of the system then the wave propagating outward will dissipate 
less than the one which will move toward the inner matter. The results indicate a rapid dissipation of the 
peak. The peak has been reduced by almost 70% at a time r = 4 fm/c. The displacement of the peak is 
determined by the sound horizon i.e. the distance travelled by the sound wave: dTCs{T)dT. We have 

taken the sound velocity, Cg = I/'s/s, independent of t for the expanding QGP background. The spatial 
anisotropic structure of the system formed in HIGRE can be understood with the help of Fourier analysis in 
terms of its various coefficients. Work on the space-time evolution of the angular power spectrum for more 
realistic initial condition for the hydrodynamical solution derived from Glauber Monte-Carlo techniques is 
under progress pS] . 

In Figs. E-E the evolutions of the spatial anisotropy with different initial geometry i.e. triangular, 
quadrangular and pentagonal for n = 3,4 and 5 respectively have been depicted. It is observed that the 
spatial anisotropies dissipate fast. The propagation of these anisotropies are primarily dictated by the velocity 
of sound in the QGP background. The propagating waves for the perturbation take shape analogous to water 
waves created on the calm surface if perturbed initially with similar geometric shape. Theoretical analysis 
of the angular power spectrum of the anisotropies will shed light on the signatures of the early stage of the 
evolving matter. 
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FIG. 2: Evolution of the spatial anisotropy with initial elliptic geometry at time tq = 0.6 fm/c (upper panel). The 
lower panel shows the geometry after a time 4 fm/c has elapsed. Hydrodynamic expansion of the QGP background 
has been taken into account. 


V. SUMMARY AND DISCUSSION 

In the present work we have studied the time evolution of spatial anisotropies for a system formed in 
HICRE. The background of the spatial fluctuations has been assumed as a thermalized expanding QGP. 
The expansion of the background has been dealt with the (2+1) dimensional relativistic hydrodynamical 
model. The anisotropy has been treated within the framework of BTE in the QGP background. The 
evolution of initial spatial anisotropies with different geometry have been studied and it is found that the 
spatial anisotropies dissipate during its propagation. Theoretical analysis of these anisotropies will help in 
understanding the early stage of the matter. We have provided an explicit relations between the fluctuations 
and transport coefficients. Angular power spectrum of the anisotropies for non-expanding background has 
been estimated in the appendix. A detail analysis of this power spectrum with hydrodynamically expanding 
background is underway. This will provide detail information on the initial anisotropies. 
Acknowledgement: G. S. acknowledges the support from Department of Atomic Energy, Govt, of India 
for this work. We would like to thank Mr Sushant Kumar Singh for providing inputs from relativistic 
hydrodynamical code. 


10 



FIG. 3: Same as Fig. but the spatial anisotropy has a triangular geometry at the initial time tq = 0.6 fm/c (upper 
panel). The lower panel shows the geometry after a time 4 fm/c has elapsed. 


VI. APPENDIX: CORRELATIONS 

In this appendix, we evaluate the correlation, (A(rii)A(n 2 )). The fluctuations, A(a:,n,t) can be written 
as: 


A{x,n,t) = 

l,m 


(42) 


with o/m = (-i)'47r/ d^nYi*^{n)A{x,n,t) and 
correlations in the following way, 


Ci{x,t)5iv5mm' Using Eq. 


42 


we can find the 


(A(f,ni,f)A(f,n2,t)) = ^ (aim-a/,^,)y'/m(ni)y;/^,(n2), 

,m' 

= X! Cl{x,t)dwSmm'Ylmini)Yif^,{h2), 

l,mj' ,m' 

^'^Ci{x,t)Yim{ni)Yi*^{n2), 

l,m 


(43) 


giving {A{x,hi,t)A{x,n 2 ,t)) = _ + l)Ci{x,t)Pi{ni ■ 712 ), since, P/(ni • 77 . 2 ) = 

defines the correlations of fluctuations observed from two differ- 
C/S. C/s are the angular power spectrum which contains the 
information of the anisotropies. Work is under progress to estimate these coefficients by solving the 


^T,n,=-iYimini)Yi*^{n 2 ). Eq. ^ 
ent directions in terms of co-efficients 


43 
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FIG. 4: Same as Fig. but the spatial anisotropy has a quadrangular geometry at the initial time tq — 0.6 fm/c 
(upper panel). The lower panel shows the geometry after a time 4 fm/c has elapsed. 


hydrodynamical equations with the initial conditions taken from Glauber Monte-Carlo method. Similarly, 
one can define these co-efficients corresponding to fc-space presentation of fluctuations. Time evolutio n of 
these co-efficients can be obtained from evolution of energy density fluctuation, A(fc,n,t) given by Eq. 
Therefore, we have 


31 


{A{x,ni,t)A{x,n 2 ,t)) = 


cPk <fk' 

(27r)3 (27r)3 


(A(fc,ni,t)A(P,n 2 ,t)) 


(44) 


Now, 


(A(fc, hi, t)A(fc', n2, t)) =(A(f, ni, to)A(T, n2, to))L{k, ffi; A/, to) 

+ (0(f, ni, to)0(T, n 2 , to))M{k, to; fc', ffi) 
-I- (A(f, ni, to)0(^, n2,to))N{k, to; k', to), 


( 45 ) 
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FIG. 5: Same as Fig.j^but the spatial anisotropy has a pentagonal geometry at the initial time tq = 0.6 fm/c (npper 
panel). The lower panel shows the geometry after a time 4 fm/c has elapsed. 


where, 




M(k, t, fc', to) =e 


k{t-to) k'{t-to) 

2 (t-tn) /40 T] smk{t — to) 3 cos k{t —to) 3sink{t — to) 


3 sT) ^ k{t-to) ' (fc(i-to))^ (Ht-to))^ 


} 


sin k'{t — to) 3 cos k'{t —to) 3 sin k'{t —to) 

k'{t-to) {k'{t-toW {k'{t-to))^ 


(46) 


2 (t-tn) 40 , 77 , .sinfcCf — in). 

N(k.t.k'M=e- . 


|.sinfc'(i — io) , 3 cos k'{t —to) 3 sin k'{t — to) ^ 

"" ^ fc'(i-io) ^ (k'it-to))^ (fc'(i-io))3 ^ 

For two functions A and 0, defining the correlation as: (A(fc,ni,i)0(fc', 77 , 2 ,i)) = {2TT)^5{k — 

A/)^Ae(A(fc,ni,i)0(fc,n2,t)) and A{k,h,t) = Y)i,m^tiik,t)Yimin), {a^im-a*f^,) = Cf^^{k,t)Su>Smm', we 
get 


21 + 1 


(A(fc, 77,1, i)A(fc', 77,2, i)) = (27r)^i5(fc - Pi{ni ■ 772 ) 


X {Cf"^{k,to)L{k,t,k,to) + Cf^^{k,to)M{k,t,k,to)}. 


00 / 


( 47 ) 
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Using {A(k,ni,t)A{k',n 2 ,t)) = 5{k- k’)'^^ -^—Cf^^{k,t)Pi{ni ■ n 2 ),in Eq.we get, 

i 

Cf^^ik, t) = C^^{k, to)L{k, t, k, to) + Cf^ik, to)M{k, t, k, to). (48) 

Using Eq. in Eq we get 

C^^{x,t)= J d^k{C^^ik,to)L{k,t,k,to)+Cf^ik,to)M{k,t,k,to)} (49) 

This provides the correlations at different angular scales at time t for a given correlations at initial time, 
to {t > to). 
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